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I. INTRODUCTION 


The properties of relativistic quantum ideal gases are a basic tool for studying more complicated equations of 
states (EoS), which very often include a sum of ideal gases. Usually, various thermodynamic quantities are expressed 
in an integral form or as a series of modified Bessel functions (see below). Unfortunately, the integrals cannot be 
evaluated exactly while the series is fast converged only in the low-temperature limit. So one has to emplo y n umerical 
integration which is quite a slow procedure. An attempt to hasten numerical calculations was made in [l|, but the 
proposed method is hardly scaled and uses the nonlinear equation solving, which is not simple either. Also, numerical 
schemes are not allowed to make qualitative conclusions while some authors need high-temperature expansion in a 
foreseeable form to explore, for example, the chiral transition in QCD 0. 

Many authors have tried to obtain the high-temperature expansion (see @,i and reference therein). The first 
terms can be easily derived (see for example [5|). For an ideal Bose gas, the task was completely solved by Haber and 
Weldon Q. Also, for the sum of a particle and its antiparticle, the high-temperature expansion of the thermodynamic 
potential, which is just the pressure up to a factor, was recently obtained in terms of special functions of a complex 
argument by Klajn [^. However, sometimes one needs to separate particles and antiparticles and so it is necessary 
to know odd power terms in the series. Also the final result for fermions calls for calculations of complex special 
functions which is not easy. So it is convenient to exclude evidently complex numbers and to give the high-temperature 
expansion in a similar form as in Q. Comparing Q and Q we see that the method used in the first paper is more 
complicated. Here we will provide how to get the same result for an ideal Bose gas and the expansion for Fermi gases 
in a simple way like in Q without complex numbers in the final expressions. 

The paper is organized as follows. In Sec. H, we introduce the integral to be studied, note how it gives the 
low-temperature expansion, and obtain the hight-temperature expansion in terms of polylogarithms, Lis(z). Our 
technique is to use the expansion of modified Bessel functions and to perform resummation of the obtained double 
series. In Sec. HI, the final result is formulated. In Sec. IV, we give the expansion for first derivatives of the pressure. 
We have included some relevant mathematical information in Appendixes A and B. Appendix C contains the proof 
of equivalence for my and Klajn’s Q results for fermions. 


II. LOW- AND HIGH-TEMPERATURE EXPANSION IN TERMS OF SPECIAL FUNCTIONS 


We consider the problem of calculation of the pressure of an ideal Bose or Fermi gas. All other quantities can be 
obtained from the pressure (see Sec. IV). 

We start with the integral representation 


P(T, /i, m) = —a T 


f >" 


/ ( 2^)3 



( 1 ) 


where E = + m? and m is the mass of particles, /x is the chemical potential, T is the temperature, jS = l/T, and 

the statistics are 


f — 1 for fermions, 
-1-1 for bosons. 


Integration of Eq. (ED by parts results in 
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It is convenient to introduce dimensionless variables: 
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Then Eq. (II} can be rewritten as 


/p(A, v) = 




2^4 


1 


LXJ 

[ dx(x2-A2)3/2 - L 


Expanding the occupation number for :/ < A so that e'^ ^ < 1, 
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Equation ([61) leads to the well-known result, 
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where Kn{x) is the modified Bessel function. This expansion is valid for any A > 0 and v < X, but it is quickly 
converging only for A ^ 1. So it is useful to have a series which can be applied at small A. 

Eor this reason we substitute the series representation of the modified Bessel function Q , 
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and using the polylogarithm function definition m with Eq. (IB3I) from Appendix B, one easily obtains the following 
common formulae, 
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which is proved for p < 0 when loe"! < 1 and is extended for any sign of p by analytic continuation. One should note 
that the expansion (nni) as a whole is not analytic. 

A private case is the expression for A = 0, 

/p( 0 ,p) =-^04(06"). ( 11 ) 

TT^ 

We also note that the nonrelativistic (low-temperature) limit Q is just a consequence of Eq. ([S]). If one uses the 
asymptotic series, 
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^ Compare with Eq. (22) from Q. 
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the result is 
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where t' = v—X is the nonrelativistic chemical potential. Using Eqs. dHU and (IB12I) . one can obtain the corresponding 
high- (e*^ <C 1) and low-temperature (e*^ ^1) expansioid^- 


III. THE HIGH-TEMPERATURE EXPANSION THROUGH ELEMENTARY FUNCTIONS. 


Equation pop immediately gives us the leading terms of the high-temperature expansion up to Eor higher 

power terms we need to deal with ^Lis(ae'^), which is not clear. So it is convenient to replace all special functions 
in Eq. (HQl) by the corresponding series. 

Eollowing [3[, we break up /p into pieces that are even and odd in p: 



/p(A, 12) + /p(A,-12) 

ip{a^) = l 

Ip{X,v) - Ip{X,-v) 

ip{x,u) = Pp{x,a + rp{x,u). 


(14) 


A. High-temperature fermions 


Then, for a = — 1, substituting Eqs. 
we get 


(IBllll . and (IB7I1 from Appendix B and changing the snmmation order, 
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^ For fermions, l/ € (— 00 , 00 ) and € [0, 00 ). For bosons, since we consider the termodynamical limit, V —>■ 00 , and do not consider the 
Bose condensation, u G (—oo,0] or G [0,1]. So there is no low-temperature expansion for bosons. 


4 





























and the hypergeometric functions F (a, b; c; are just polynomials: 
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F (a, b; c; are polynomials since a and b are negative integers or zero. 

It can be shown that Eq. m completely coincides with the result Q for fermions. The corresponding proof is 
given in Appendix C. 

Rewriting F (—fc, —k — 2; i; r^) as Jacobi polynomials and applying the recurrence relation, we find that the series 
representation of Ip{X, v) is converged for A + |i/| < tt. 


B. High-temperature bosons 


In the same way, for a = 1, using Eqs. 
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which exactly coincide with the result (3|. Here, 


b{x) = r(a:)C(x). 


This exercise allows us to be sure of the correctness of the result for the Fermi-Dirac integral. 
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IV. OTHER THERMODYNAMICAL QUANTITIES 


Besides the pressure, three other quantities are often used: the particle number density, the scalar density, and the 
entropy density which are the first derivatives of the pressure. 
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s{T,fi,m) = [4-fp(A,;/) + XlsdX^v) - p/„(A, i/)] 

= r3/,(A,p). 

As a result, from the first law of thermodynamics, one obtains the energy density 

e(r, /i, m) = Ts{T, /i, m) + /i n{T, ji, m) — P(T, /x, m) 
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Below we list the complete high-temperature expansion of the Fermi-Dirac integrals for the particle number density, 
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All functions G^(r) are polynomials. 


V. SUMMARY 

The high-temperature (A —?► 0) expansion of the integral ([T]) for Fermi-Dirac statistics has been obtained. For 
thermodynamic applications, this result allows us to obtain all thermodynamic quantities. The proposed method for 
obtaining high-temperature expansion is simple and reproduces the results of [3|, 1^ . 
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Appendix A: Properties of gamma, digamma, and Riemann zeta functions 

In our consideration, we often meet three special functions: gamma, digamma, and Riemann zeta functions. In this 
Appendix we just quote some useful relations for these functions. 

First of all, we have @ 
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Then applying Eq. (IA2I) to the nominator and denominator, we obtain 
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Also we permanently keep in mind that 


r(n) = (n-l)!, r(l/2) = V^, 
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where n G N. 

In our computations we need to know how to take the derivative of the gamma function. So we deal with 

r'(x) = '>p{x)r{x) 

and we use the following features of the digamma function Q , 
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where je is Euler’s constant. From Eqs. (ITtI) and (IA2I) it can be found that 
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The last remarkable object which occurs in the paper is the Riemann zeta function, (^{x). We use that @ 

C(z) = 2(27r)^“^ r(l - z) C(1 - z) sin ^ 


(A6) 


(A7) 

(AS) 


(A9) 


(AlO) 


and 

C(0)=-i, C(-2n)=0, = 

C'(0) = ln(27r), 

where n G N. 
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Also differentiating both sides of Eq. (lAlOp . we get 
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for n G N. 

For brevity, we sometimes also use the Dirichlet eta function 

r^iz) = (1 - 2^-^) C(^) (A15) 

whose properties follow from the properties of C(2^)- 
























Appendix B: Properties of polylogarithms 


The polylogarithm function can be defined by power series Q, 
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We need a series expansion of Lis(e^) and Lis(—e^) for small z. The first one is obtained in 3, and we quote the 
result: 
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where Eq. (lAllI) is used. 

To obtain the expansion for Lis(— e^), one should just apply Eq. (IB4I) for z + itt or use Taylor’s theorem and Eq. 
|. Then it results in Q 
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For convenience, we also separately give the expressions for s = —2m, m = 0,1, 2,.. 


/ zx r( 2 m + 1 ) 1 

Fl-2m(e j — ^2m+l 2 


( 27 r) 2 m+l 


E (-1) 


^T{2k)Ci2k) / Z \2k-2m-l 






Li — 2 m ( ^ ) — 2 ‘^'^0 


-2m+l 


E (-1) 


r(2fc- 2m) 

fcr(2A:)C(2fc)(l - 2-2'=) ,z\ 2 k- 2 m-l 


k—Tn-\-l 


r(2fc - 2m) 




(B8) 


(B9) 


which are directly obtained from Eqs. (IB4I) and (IB6I) where Eq. (lAllD is taken into account. 

Finally, we should find the derivatives of polylogarithms over the index for s = —2m. It can be easily made by 
taking the corresponding derivative on the rhs of Eq. (lB4l) . 
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To obtain the last two expressions, one should use Eqs. (IA13I) and (IA14E 
Let us also note the asymptotic expansion 
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Appendix C: EQUIVALENCE OF EQ. ([15]) TO KLAJN’S RESULT 0. 

Applying Eq. (51) of 0 for fermions so that Pfl = v + in, we obtain 
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We need to prove that this expression is equal to Eq. (unj. 

The Bernoulli polynomials are Bi{z) = z'^ — 2z^ + — 1/30 and B 2 {z) = — z + 1/6. This gives 
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The derivatives of the polygamma functions may be expanded in a Taylor series about the argument 1/2: 


,<«. (i_ -)(i, -) = (i) 




(C4) 


Collecting the last three expressions, separating out the term 2k + 2j = 0, and changing the sum order, results in 
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for I G N. 

The last formulas together with 1 / = rX immediately lead Eq. (jC5l) to the form of Eq. (Ha. 
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